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Abstract 

The cross sections for exclusive muon pair production in nucleus - nucleus collisions are calculated 
and several differential distributions are shown. Realistic (Fourier transform of charge density) 
charge form factors of nuclei are used and the corresponding results are compared with the cross 
sections calculated with monopole form factor often used in the literature and discussed recently 
in the context of higher-order QED corrections. Absorption effects are discussed and quantified. 
The cross sections obtained with realistic form factors are significantly smaller than those obtained 
with the monopole form factor. The effect is bigger for large muon rapidities and/or large muon 
transverse momenta. The predictions for the STAR and PHENIX collaboration measurements at 
RHIC as well as the ALICE and CMS collaborations at LHC are presented. 
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I. INTRODUCTION 



In Fig{T]we show the basic QED mechanism of the exclusive production of muon pairs. 
The shaded circles represent the coupling of photons to large-size objects - nuclei. In the 
momentum space this is done in terms of electromagnetic form factors of nuclei. In the case 
of scalar nuclei there is only one form factor - the charge form factor of the nucleus. 



A2 




FIG. 1: (Color online) The Born diagram for the exclusive dimuon production. 

It was recognized long ago that the production rate of leptons in ultrarelativistic heavy 
ion collisions is enhanced considerably by the coherent effects and large charge of colliding 
ions [1]. Many results have been presented in the literature since then (for reviews of the 
field see e.g. jH, Isj). Recently, there was a growing theoretical interest in estimating higher- 
order QED corrections j#i7|. In most of the practical calculations of exclusive dilepton 
production a simple monopole charge form factor of the nucleus was used. While it may be 
sufficient for estimating the total cross section, it may be not sufficient for calculations of 
the differential cross sections. The importance of including realistic charge form factors was 
discussed recently for exclusive production of pairs of p° mesons j8|. 

Most of the existing calculations concentrated on total cross section, interesting theoreti- 
cal quantity, which cannot be, however, measured in practice, neither at RHIC nor at LHC 
The experiments running at RHIC and those planned at LHC demand severe cuts on lepton 
transverse momenta or on their rapidities. 

It is the aim of the present analysis to make realistic estimates of the cross sections 
including the experimental cuts. We shall compare the results obtained with monopole form 
factor used in the literature and the results obtained with realistic form factor being Fourier 
transform of the charge density of the nucleus. We shall perform the calculation in the 
equivalent photon approximation (EPA) in the impact parameter space as well as in the 
momentum space. While the impact parameter EPA allows to include easily absorption 
effects due to the size of colliding nuclei, the momentum space approach allows to study 
easily several differential distributions. In our calculation we shall include experimental 
limitations of the STAR and PHENIX detectors at RHIC and those of the ALICE and CMS 
detectors at LHC. 



2 



II. FORMALISM 



A. Charge form factor of nuclei 

The charge distribution in nuclei is usually obtained from elastic scattering of elec- 
trons from nuclei jof. The charge distribution obtained from these experiments is often 
parametrized with the help of two-parameter Fermi model (lol |: 

p(r) = po (l + exp(^^)) , (2.1) 

where c is the radius of the nucleus and a is the so-called diffiusness parameter of the charge 
density. 
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FIG. 2: (Color online) The ratio of the charge distibution (p) to the density in the center of nucleus 
iPo). 

Fig. [2] shows the charge density normalizationed to unity at r = 0. The correct normal- 
ization is: pAn(O) = ^-^^^ fm~^ for Au nucleus and ppb{0) = °'^^°^ /m~^ for Pb nucleus. 

The form factor (F) is the Fourier transform of the charge distribution [9| . If p (r) is 
spherically symmetric then the form factor is a function of photon virtuality (g) only: 

F(q) = / — p (r) stn (qr) rdr = 1 - + .... (2.2) 

J q 3! 5! 

Fig. [3] shows the moduli of the form factor as a function of momentum transfer. The 
results are depicted for the gold (solid line) and lead (dashed line) nuclei for realistic 
charge distribution. The realistic form factor is obtained as a Fourier transform of the 
realistic charge density which we take from the literature jof. Here one can see many 
oscillations characteristic for relatively sharp edge of the nucleus. For comparison we show 
the monopole form factor often used in the literature. The two form factors coincide only 
in a very limited range of q and with larger value of q the difference between them becomes 
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FIG. 3: (Color online) The moduli of the 
charge form factor {q) of the ^^"^Au and 
208p^ nuclei for realistic charge distributions. 
For comparison we show the monopole form 
factor for the same nuclei. 
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FIG. 4: (Color online) The monopole form 
factor for the values of A reproducing charge 
radius of ^^'^Au and ^osp^ nuclei and for com- 
parison for A = 0.08 GeV often used in the 
literature. 



larger and larger. 



The monopole form factor [HI] given by the simple formula: 



(2.3) 



A2 + g2 

leads to a simplification of many formulae for photon-photon collisions. In our calculation 



A is adjusted to reproduce the root mean square radius of a nucleus (A 



with the 



help of experimental data [10|: 

• for ^^'^Au: < >^/'^= 5.3 ^ A = 0.091 GeV, 

• for P6: < r2 >i/2= 5.5016 ^ A = 0.088 GeV. 

Different values of A are used in the literature, ranging from 80 to 90 MeV. Fig. H] shows the 
monopole form factor with A adjusted to reproduce the rms radius of the charge distribution. 



B. Equivalent Photon Approximation 

The equivalent photon approximation is the standard semi-classical alternative to the 
Feynman rules for calculating cross sections of electromagnetic interactions 12j. This is 
illustrated in Fig. |5] where we can see a fast moving nucleus with the charge Ze. Due to 
the coherent action of all the protons in the nucleus, the electromagnetic field surrounding 
(the dashed lines are lines of electric force for a particles in motion) the ions is very strong. 
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FIG. 5: (Color online) Equiva- FIG. 6: (Color online) The quantities 
lent photon approximation. used in the impact parameter calcula- 

tion. 

This field can be viewed as a cloud of virtual pliotons. Tliese pliotons are often considered 
as real. They are called "equivalent" or "quasireal photons". In the collision of two ions, 
these quasireal photons can collide with each other or with the other nucleus. So the strong 
electromagnetic field is used as a source of photons to induce electromagnetic reactions on 
the second ion. We consider very peripheral collisions. It means that the distance between 
nuclei is bigger than the sum of the radii of the two nuclei {b > Ri + R2 = 14/m). Fig. [6] 
explains the quantities used in the impact parameter calculation. We can see a view in the 
plane perpendicular to the direction of motion of the two ions. In order to calculate the 
cross section of a process it is convenient to introduce a new kinematic variable: x = J^, 
where uj is the energy of the photon and the energy of the nucleus Ea = 'jArriproton = 7^a, 
where is the mass of the nucleus and 7 is the Lorentz factor. 
The total cross section can be calculated by the convolution: 

a (^AA n'^ji'AA; saa) = J ^ (^7 A*^/^"; ^-y-y = V^i^^saa) dn^j (xi, X2, b) . (2.4) 

The effective photon fluxes can be expressed through the electric fields generated by the 
nuclei: 

1 1 
dra^^ (a;i,a;2,b) = -d^bi|E (xi, bi) p-d^b2|E (0:2, b2) P 



TT TT 



X 



^L(b)5(^)(b-b, + b2)^^. (2.5) 



The presence of the absorption factor S*^^^ (b) assures that we consider only peripheral 
collisions, when the nuclei do not undergo nuclear breakup. In the first approximation this 
can be expressed as: 

51 (b) = 9{h- 2Ra) = 9 (|bi - b2| - 2Ra) . (2.6) 

Thus in the present case, we concentrate on processes with final nuclei in the ground state. 
The electric field strength can be expressed through the charge form factor of the nucleus: 

E (X, b) = ZV^^ J (0) ^'^' ^2 +^2m/ - + ^'^^) • ^2.7) 
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Next we can benefit from tfie following formal substitution: 

i J d%\E (a;, b) |2 = J d^N {to, h) = n (cu) (2.8) 

by introducting effective photon fluxes which depend on energy of the quasireal photon u 
and the distance from the nucleus in the plane perpendicular to the nucleus motion b . 
Then, the luminosity function can be expressed in term of the photon flux factors attributed 
to each of the nuclei 

dn^^{ui,U2,h) = J 9{\hi-h2\-2RA)N {uji,hi)N {uj2,h2)(i'^^id'^^2dujidu2. (2.9) 

The total cross section for the AA — t- fi^fi^AA process can be factorized into the equivalent 
photons spectra ( n {u) ) and the 77 subprocess cross section as (see e.g.[13|): 

a (^AA saa) = J cr (^-fj l^^l^'] ^"77) 6* (|bi - ba] - 2Ra) 

X (wi, bi) {002, ba) d^bid^bgdwidwa , (2.10) 

where W^^ = a/IcJio^ is energy in the 77 subsystem. Eq. ( I2.10p is a generalization of the 
simple parton model formula (see e.g. [2|): 



a (^AA — )■ /i^/i AA^ ~ J ^ (^^ ~^ ' n i^i) n (1^2) dwidcja • (2-11) 

Additionally, we define ^ = | (y^+ + y^- ) , rapidity of the outgoing dimuon system which is 
produced in the photon-photon collision. Performing the following transformations: 

c^i = 0:2 = ^e-^ , (2.12) 

dwidwa = ^dW^^dr , (2.13) 



dwidwa — ^ dPV^^dF where 
formula f l2.10p can be rewritten as: 



d{ui,U2) 



d{W,,,Y) 



(2.14) 



a [AA ij+fi~AA; saa) = y a (77 ^ f^^l^ ; ^^77) ^ (|bi - bsl - 2Ra) 



xN{uJi,hi)N{uj2,h2)^d%id^h2dW^^dY . (2.15) 
Finally, the cross section can be expressed as the five-fold integral: 

a (^AA n+fi-AA; saa) = J ^ {ll ^ l^^l^'; W^^) (|bi - bsl - 2Ra) 

xN {uJi,hi)N {u2,h2)27ib^db^dKdby^^dW^^dY , (2.16) 

where bx = (&ix + &2a;)/2, by = {biy + b2y)/2 and bm = bi — 62 have been introduced. This 
formula is used to calculate the total cross section for the AA — )■ AAfi~^fi~ reaction as well 
as the distributions in 6 = bm, W^^ = M^+^- and F(/i+/x~). 
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Different forms of form factors are used in the literature. We compare the equivalent 
photon spectra for an extended charge distribution (realistic case) to the monopole case. 
The dependence of the photon flux on the charge form factors can be found in 0|: 



N {co, h) 



I 



62 



J nVi (u) F 



( 



V 



V 7 



+ m2 



62 



V 



-du 



(2.17) 



where J\ is the Bessel function of the first kind and q is the four-momentum of the quasireal 
photon. The calculations with the help of realistic form factor are rather laborious, so often 
a simpler monopole form factor is used [lH. Introducing monopole form factor to f l2.17p one 
gets: 



A^(w,6) 



^ "em 1 



bu 
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(2.18) 



where Ki is the modified Bessel function of the second kind. 




FIG. 7: (Color online) The equivalent photon number as a function of impact parameter (integrated 
over w), see Eq. (j2.19p . 



Fig. [7] shows the distribution of the equivalent photon number as a function of the impact 
parameter 

N{b) = J N{u, b)du. (2.19) 

We present the results for gold and lead nuclei, for realistic and monopole form factors. Here 
we do not impose any sharp cutoff on the impact parameter. One can see that for small b 
the flux factor with monopole form factor is bigger. For large b the results obtained with 
the help of realistic and monopole form factors are almost the same. 

In addition in Fig. [8] we show the ratio of equivalent photon fluxes obtained with the help 
of realistic form factor to that for the monopole form factor. The oscillations in b are due to 
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FIG. 8: (Color online) The ratio of the flux factor obtained with realistic charge distribution to 
that with the monopole form factor as a function of impact parameter. 
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FIG. 9: (Color online) The equivalent photon number n (w), see Eq. ()2.20p . Left panel: h G (0, 100) 
fm, right panel: 6 G (14, 100) fm. 



step-like distribution of the charge in the nucleus. The results for lower (v^atat = 200 GeV) 
and higher (a/sjva? = 5.5 TeV) energies are almost the same. 
Fig. |9] shows 

N{uj)= [ 2'KhN (w, h) db. (2.20) 



Here we consider the integral over full range of the impact parameter (left panel) and for 
b > 2Ra (right panel). One can see that the difference between monopole and realistic form 
factor for both gold and lead nuclei is not significant. The quantity shown depends rather 
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weakly on the photon energy. 
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FIG. 10: (Color online) The elementary cross section for the 77 — t- ^^/U reaction as a function of 
the photon-photon energy. 

In Fig. [To] we show the energy dependence of the elementary 77 — )■ yU^/i^ cross section 
used in our EPA calculations [l|]: 



a (77 ^ ) 



77 

X J 2 In 



(2.21) 



77 



2m„ 



'l + v) 



K ^77 



^ ^77 ' 



where 



^ 77 

This formula is often called the Breit- Wheeler formula. 



Ami 

1 _ 



(2.22) 



C. Momentum space calculation 

We consider a genuine 2 — )■ 4 reaction (see Fig. [TT]) with four-momenta Pa + Pb 
Pi + P2 + P3 +Pa- In the momentum space approach the cross section for the production of 
a pair of particles can be written as: 



cr 



— I |2 {2Tlf 5^ {pa +Pb-Pl-P2-P3- Pi) 



^ d^Pi d^p2 d^Ps d^PA .2 23) 

{2Tif 2Ei (27r)^ 2E2 (27r)^ 2E3 (27r)^ 2E4 ' 
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FIG. 11: (Color online) Amplitude of the considered process. On the left one can see the t-channel 
amplitude and on the right - the u-channel amplitude. 



Using 



dyi d^pit = dyiPitdp, 



(2.24) 



Eq. fl2.23p can be rewritten as: 

^ = J ^\M\^^'^{Pa+Pb-Pl-P2-P3-Pi) ,8 



1 1 



X 



(27r)' 24 

{dyipitdpitd(j)i) {dy2P2tdp2td(j)2) (dysd'^pst) {dyd^pii^ . 



(2.25) 



In the above formula pu are transverse momenta of outgoing nuclei and leptons, (pi, 02 are 
azimuthal angles of outgoing nuclei. Additionally, we introduce a new auxiliary quantity 

Pm = P3i - P4t (2.26) 

and benefitting from 4-dimensional Dirac delta function properties, Eq. f l2.25p can be written 
as: 



cr = J yJM\^ (Ea + Eb- El- E2- Es- E^) 5^ {pu + P2z + P3z + Paz 



1 1 



(27r)' 24 



X {dyipitdpitd(f)i) {dy2P2tdp2td(t)2) dy^dy^d^Pm ■ 



(2.27) 



The energy-momentum conservation gives the following system of equations that has to be 
solved for discrete solutions 



/i - ^3 - ^4 = 
-PZz - PAz = Plz + P2z , 



mlt+piz + ^mlt+plz , 



(2.2^ 



where mu, m2t are the so-called transverse masses of outgoing nuclei which are defined as: 

ml=pl + ml. (2.29) 

We wish to make the transformation from (?/i, 1/2) to {piz-, P2z)- The transformation 
Jacobian takes the form: 



jTfc 



Plz {k) 



P2z (k) 



rriit + Pizik) Jmlt+pl,{k) 



(2.30) 
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where k numerates discrete solutions of Eq. (12.281) . Thus the cross section for the 2 — )■ 4 
reaction reads: 

1 ^11 



= fY.^k^ (pit, 01, P2t, 02, 1/3, 1/4, Pm, 0m) , \M\ 



2^s {s - 4m2) (27r)'24 

X {pudpitd(f)i) {P2tdp2td(j)2) ^dy^dy^d^pm ■ (2.31) 

For photon-exchanges, considered here, it is convenient to change the variables pu = 
logio (Pit), P2t ^2 = logio iP2t)- The lepton helicity dependent amplitudes of the process 
shown in Fig. [11] can be written as: 

A^A3,A, (t-channel) = eF^^Qi) (Pa + PiT ^^u{ps, A3) ^YlM:L^±^ 

X lYv (p4, A4) {Pb + P2f e F,^ (gs) (2.32) 

qi + IE 



and 



A^A3,A4 (^-channel) = eFch{qi) {Pa + PiT ^^^^u{P3, A3) ^^^^ ^) + "^a*] 



2 I • \Foi ^M/ / / \2 o 

gf + «£ fe-Ps) -"^^ 

X 2 t; (p4, A4) (P6 + P2)^ e F,, (ga) • (2.33) 

'?2 

These amplitudes are calculated numerically. Finally, to calulate the total cross section 
one has to calculate the 8-dimensional integral inserting J^x^^x^^ = M A3,A4 (t-channel) + 
A^A3,A4 (u-channel) into Eq. (12.311) . We shall compare the impact parameter EPA results 
with the exact ^ Quantum Electrodynamics results. 



III. RESULTS 

Let us start from the presentation of the results obtained in the impact parameter EPA. 
In FigJT2] we show the distribution in the impact parameter b for typical RHIC energy 
^/sNN = 200 GeV. The contributions from distances smaller than b = 2Ra are cut off and 
consistently with ^-function in Eq. (I2.16p . We clearly see a huge contribution from distances 
large compared to the nuclear size. The distribution with realistic charge falls off somewhat 
quicker than that for the monopole charge form factor. This is better visualized in the right 
panel where the ratio of the corresponding cross sections is shown. 

The difference of the cross sections for the monopole and exact charge form factors at 
large impact parameter b shown in the figure is especially intriguing in the light of the 
equality of the photon flux factors at large bi or 62 (see FigJT]). How to understand this 
quite nonintuitive result? In FigJiniwe show the distribution of da/dbidb2 in (&i,&2) with 
the severe restriction for the impact parameter b G (480,520) fm. We see two pronounced 
peaks at (61 ^ 6, 62 ~ 0) and (61 ^ 0,62 ~ b). This demonstrates a strong preference of 



By exact we mean the correct inclusion of the 2 — )■ 4 process phase-space. It is, however, rather difficult 
to include absorption effects in this approach. 
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FIG. 12: (Color online) The cross section as a function of the impact parameter for the AuAu — >■ 
AuAu reaction calculated in the equivalent photon approximation. In the left panel we show 
the results for realistic charge distribution (solid line) and for monopole form factor (dashed line). 
On the right side we depict the ratio : RATIO = da (f^^^^^^^^^^ /da {P^i^^OPOLE^^ _ 




asymmetric production of the pair: close to the trajectory of one or the other nucleus, where 
the form factor details are important (see FigJTj). This point was never discussed so far in 
the literature. 

The distributions shown in Fig{T2] are purely theoretical, that is cannot be easily mea- 
sured. Let us come now to the distributions which could, at least in principle, be mea- 
sured. Fig{T4] shows the distribution in the dimuon subsystem energy. The distributions in 
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W^-y = Mfj_+fj_- falls steeply off. In the right panel we show the ratio of the cross sections for 
realistic charge distribution to that for the monopole charge form factor. At W-y-y = 10 GeV 
the two distributions differ already by a factor of about 5 which clearly shows limitations of 
the calculations with analytic charge form factors. 




FIG. 14: (Color online) The cross section for Au — Au scattering as a function of photon-photon 
center-of-mass energy Wyy = M^+^- in EPA. In the right panel we show the ratio of "realistic" 
to "monopole" form factor. 

Finally, in analogy to the AA — )■ AAp^p^ reaction studied in Ref. [s!], in FigJl5]we show the 
distribution in the dimuon pair rapidity. As for the p°p° production we see a huge difference 
between the results of the two calculations for large dilepton rapidities. Measurements of 
dileptons in forward directions would be therefore very useful to understand the role of 
realistic charge distribution. The relative effect is shown in the right panel of the figure. 

The preliminary calculation in the impact parameter space clearly shows how important 
can be studying of differential distributions to pin down the effects of realistic charge density. 
Not all of the distributions can be easily addressed in the impact parameter approach. The 
Feynman diagram approach in the momentum space seems to be a better alternative to 
study the differential distributions. 

Now we come to the presentation of results obtained in the momentum space approach 
with details outlined in Section II. FigJTGl shows distributions in muon rapidities (identical 
for /i+ and yU~). No other limitations or kinematical cuts have been included here. As in the 
previous cases we show distributions obtained with the monopole and realistic charge form 
factor. The effect of the oscillatory character of Fch{q) and in particular its first minimum 
is reflected by a smaller cross section at larger rapidities compared to the results obtained 
with monopole form factor. This is due to the fact that on average at large rapidities larger 
four- momentum squared transfers [ti or ^2) are involved. In reality, one effectively integrates 
over a certain range of ti and ^2- The relative effect is shown in the right panel. 

Fig JT71 shows the situation (the ratio of the two calculations) in the two dimensional space: 
(7/3,1/4). Clearly at mid rapidities, where on average rather small ti and ^2 are involved, the 
use of the approximate monopole form factor is justified. This is not the case at the edges 
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FIG. 15: (Color online) The cross section as a function of y = ^ {y^+ +y^-) (left panel) for 
realistic and monopole form factors (left) calculated in EPA and their ratio (right). 
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FIG. 16: (Color online) The cross section as a function of y^- for realistic and monopole form 
factor calculated in the momentum space (left panel). Their ratio is shown in the right panel. 



of the (i/3,?/4) plane where due to kinematics |ti| or/and |t2| are larger. 

Up to now we have discussed " a theoretical situation" when all the muons are accepted. 
In practice one can measure only muons with transverse momenta larger than a certain 
value, characteristic for a given detector. We shall consider now cases relevant for concrete 
experimental situations. 

The calculations in the literature concentrated mostly on the total cross section. In 
FigiH] we present the dependence of the total cross section on the lower cut-off in the 
impact parameter. We present EPA results for realistic (lower solid line) and monopole 
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FIG. 17: The ratio of two-dimensial distributions da (f^^^^^^'^^^) /da (f^^^^^^^^^ in and 
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FIG. 18: (Color online) The compilation of the results obtained in different approaches for the 
total cross section for AuAu — t- AuAu^^ at -v/stva^ ~ ^^O GeV. 

(upper solid line) form factors. The cross section without the cut-off is by 15% larger than 
that for hcut = 14 fm. This result is smaller than the corresponding results obtained within 
momentum space calculations, shown as the horizontal dashed lines. Different methods has 
been used in the literature to calculate the total cross section for the AuAu — )■ AuAufi'^ fi~ 
process. For comparison we show also results obtained recently by Jentschura and Serbo 
(JS) [l^ in the momentum space EPA and by Baltz et al. jl5| in the b-space EPA. The JS 
result should be compared to our momentum space calculation with monopole form factor. 
Our exact calculation is in this case larger than their EPA calculation by about 24%. This 
shows the precision of the momentum space EPA. The Baltz et al. result is significantly 
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lower than our b-space EPA result. In their calculations the cuts were imposed rather on bi 
and 62, instead on b in our case. If we impose additional cuts on bi and 62 in Eq- f !2.16p we 
get the point in the lower-right corner. If the cut on b is not imposed we get the point in 
the lower-left corner. The result of Baltz et al. differs from both these values, the solution 
being most probably a different form factor used in their case. 

Now we will continue reviewing our predictions for the differential distributions. Let us 
start with the ALICE detector. The ALICE collaboration can measure only forward muons 
with psudorapidity 4 < 77 < 5 and uses a relatively low cut on muon transverse momentum, 
Pt > 2 GeV. In Fig{T9] (left panel) we show the invariant mass distribution of dimuons for 
monopole and realistic form factors. The ALICE experimental cuts were incorporated into 
our calculations. The bigger invariant mass the bigger the difference between the results for 
the two form factors. The same is true for distributions in muon transverse momenta (see 
the right panel). 




FIG. 19: (Color online) Invariant mass distribution ^^j^f^ — (left) and muon transverse momentum 

distribution ^ = ^ (right) for ALICE conditions: 2/3,2/4 G (3,4), pst^pu > 2 GeV and the 
center-of-mass energy Wnn = 5.5 TeV. 

The distribution in rapidity is shown in FigEUl We present the cross sections for both 
(realistic, monopole) form factors and their ratio. 

Double differential distribution of the muon rapidity and transverse momentum is shown 
in Figj2Tl These are our predictions which could be studied experimentally in the future. 
The small irregularities seen in the two-dimensional spectra for realistic form factor are the 
consequence of the oscillatory character of the nucleus charge form factor. The distribution 
for the monopole form factor is more smooth. 

In FigJ22] we show the ratio of the cross sections shown in the previous figure. Huge 
deviations from the unity can be seen. The reminiscence of the oscillating form factor can 
be seen also in the ratio. Experimental confirmation of this behaviour would be very useful. 
Moreover it would demonstrate whether our understanding of the nuclear effects is correct. 
Large deviations from the predictions presented here would be surprising. 

Let us come now to the predictions for the CMS detector. In contrast to the ALICE 
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FIG. 20: (Color online) ^ = ^ (left) and their ratio (right) for the ALICE conditions: 7/3,7/4 G 
(3,4), p3t,pit > 2 GeV and Wnn = 5.5 TeV. 





FIG. 21: Double differential cross section ^^^^ for realistic (left) and monopole (right) form 
factors for ALICE conditions y3,y4 G (3,4), pstiPu > 2 GeV and Wnn = 5.5 TeV. 



detector, CMS can measure midrapidity values with -2.5 < 7/3,7/4 < 2.5. At midrapidities 
one samples on average smaller ti and t2 therefore the efects of the realistic form factors 
are expected to be smaller. FigJSn] confirms the expectations. Even for muon transverse 
momenta of 50 GeV one obtains damping with respect to the result obtained with the 
monopole form factor by a factor of about two only. 

The cross section dependence on the muon rapidity is shown in Fig|24l Rather large 
cross section of the order of 0.1 mb is expected within the CMS acceptance. The average 
deviation with respect to the monopole form factor is about 20% (see the left panel). 

The two-dimensional distributions within the main CMS detector are shown in Fig{ 
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FIG. 22: Ratio of the cross sections ^1/3 dpat ALICE conditions: 2/3,1/4 G (3,4), P3t,P4t > 2 

GeV and PFjviV = 5.5 TeV. 
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FIG. 23: (Color online) The muon transverse momentum distribution (left) and the ratio 
(right) for the CMS conditions: 2/3,2/4 G (-2.5,2.5), P3t,Pit > 4 GeV and Wnn = 5.5 TeV. 



Big modifications with respect to the monopole case can be seen for large pt and \y^+,yfj,- \ ~ 
2.5, which is also presented in the form of the ratio in Figj26l 

Finally, for completeness in Figj27]we show the distributions in the {yziVi) plane. Here 
the distributions obtained with the monopole and realistic form factors are rather similar, 
but one should realize that these distributions are dominated by muons with small transverse 
momenta that are only slightly bigger than the experimental acceptance pt > 2.5 GeV and 
as a consequence relatively small ti and ^2 values. 

The same processes can be also studied at the being presently in the operation RHIC. 
Here STAR and PHENIX detectors can be used. The distribution of the muon transverse 
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FIG. 24: (Color online) The muon rapidity distribution (left) and the ratio (right) for the CMS 
conditions: 1/3,2/4 G (—2.5,2.5), P3t,P4t > 4 GeV and Wnn = 5.5 TeV. 




FIG. 25: dysdpst realistic (left) and monopole (right) form factors for the CMS conditions: 
2/3,2/4 G (-2!5.2,*5), pst,Pu > 4 GeV and Wnn = 5.5 TeV. 

momentum is shown in FigJ2Sl The STAR rapidity cuts -1 < ?/3, 2/4 < 1 are taken into 
account. Compared to the LHC the transverse momentum distributions decrease much 
faster. This fast fall-off limits the real measurements to relatively small transverse momenta 
of the order of 10 GeV. The inclusion of realistic charge distribution is here much more 
important than for the CMS conditions. The relative effect of damping with respect to the 
results with the monopole charge form factor is shown in the right panel. At pt = 10 GeV 
the damping factor is as big as 100! Experiments at RHIC have a potential to confirm this 
prediction. 

In general, one could also inspect the rapidity distributions. Our predictions are shown 
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FIG. 27: for realistic (left) and monopole (right) form factors for the CMS conditions: 

2/3,2/4 G (-2^5, 2*5), ps^Pu > 4 GeV and Wnn = 5.5 TeV. 

in Figj29l We predict the 30-40 % cross section damping with respect to the reference 
calculation (monopole charge form factor). 

The two-dimensional distributions in muon rapidity and muon transverse momenta are 
shown in FigjnU] for the realistic and monopole form factors. Their ratio is presented in 
FigJSTl Again as for the transverse momentum distribution (see Fig j28|) a huge damping 
can be observed. The irregular structure of the ratio reflects the strong nonmonotonic 
dependence of the charge form factors of Au nuclei on ti and ^2- For completeness in 
Figj32] we show the distribution of the dimuon invariant mass. The effect of the form factor 
oscillations shows up at large dimuon invariant masses where the "realistic" cross section is 
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FIG. 28: (Color online) ^ (left) and the ratio (right) for the STAR conditions: 2/3,1/4 G (-1, 1), 
P3t,Pu > 1 GeV and Wnn'= 200 GeV. 
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FIG. 29: (Color online) (left) and the ratio (right) for the STAR conditions: 2/3,2/4 
P3t,Pu > 1 GeV and Wnn = 200 GeV. 



e (-1,1), 



rather small. 

The PHENIX collaboration can measure muons in a rather limited range of rapidities 
shown in Figj33l We have given names to the four possible regions (squares) in the figure. 
In spite of these limitations, still interesting measurements can be done. As an example 
in FiglSH and Figj35]we show our predictions for SQUAREl and SQUARE2, respectively 
(the results for SQUARES and SQUARE4 are not shown as can be obtained by symmetry). 
Again large deviations from the monopole form factor results are predicted. 
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FIG. 30: dysdpst realistic (left) and monopole (right) form factor for the STAR conditions 
y3,y4 e (-M)%3t,^?4t > l GeV and Wnn = 200 GeV. 




FIG. 31: The ratio of the two-dimensional distributions from the previous figure for the STAR 
conditions: 2/3,^4 G (—1, 1), pzf,Vu ^ 1 GeV and Wnn = 200 GeV. 

IV. CONCLUSIONS 

The production of charge leptons in heavy ion collisions was proposed recently as a "lab- 
oratory" for studying Quantum Electrodynamics effects, in particular the multiple photon 
exchanges. While very interesting theoretically it is still nonrealistic because of other ap- 
proximations made in the calculations. 

In this paper we have presented a study of the role of charge density for the differential 
distributions of muons produced in exclusive ultra-peripheral production in ultrarelativistic 
heavy ion collisions. Most of the calculations in the literature use so-called monopole charge 



22 




10 20 30 40 

M^,^. (GeV) 



FIG. 32: (Color online) Invariant mass distribution g^g^^ — for the STAR conditions: 2/3,2/4 G 
(-1, 1), P3t,Pu > 1 GeV and Wnn = 200 GeV. 
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FIG. 33: (Color online) The muon rapidity regions available by the PHENIX detector. 



form factor, which allows to write several formulae analytically. While it may be reason- 
able for the total rate of the dimuon production it is certainly too crude for differential 
distributions and for the cross sections with extra cuts imposed on transverse momenta of 
muons. 

We have performed calculations in the Equivalent Photon Approximation in the impact 
parameter space and in the momentum space using Feynman diagrammatic approach. The 
first method is very convenient to include absorption effects, while the second one allows to 
study differential distributions. 

Our calculations show that the results obtained with the realistic and the approximate 
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FIG. 34: (Color online) SQUARE 1: ^ (left) and ^ (right) for the PHENIX conditions: 1.2 
< 1^3, y4| < 2.4, p3t,Pit > 2 GeV and Wnn = 200 GeV.' 




FIG. 35: (Color onhne) SQUARE 2: ^ (left) and ^ (right) for the PHENIX conditions: 1.2 
< 1^3, y4| < 2.4, p3t,Pit > 2 GeV and Wnn = 200 Gev!* 



form factors can differ considerably, in some parts of the phase space even by orders of 
magnitude. The effects related to the charge distribution in nuclei are particularly important 
at large rapidities of muons and at large transverse momenta of muons. 

We have also discussed the role of absorption effects which can be easily estimated in the 
impact parameter space. This allows to estimate the absorption effects for the total rate 
or for the rapidity distribution of the dimuon pairs. Estimating this effect in the case of 
differential distributions is not simple, but could be studied in the future. 

We have presented predictions for the STAR and PHENIX detectors at RHIC as well as 
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for the ALICE and CMS detectors at LHC. In all cases we have found significant deviations 
from the reference calculation for the monopole form factor. It would be interesting to pin 
down the effects discussed here and verify the present predictions in future studies at LHC. 
Both ALICE and CMS detectors could be used in such studies. 

In practice such studies may not be simple as an efficient trigger for the peripheral 
collisions is required. The multiphoton exchanges leading to additional excitation of nuclei 
and subsequent emission of neutrons could be useful in this context (see e.g. [Hj). The 
neutrons could be then measured by the Zero Degree Calorimeters. First measurements of 
this type for e'^e~ pair emission have been already performed by the STAR and PHENIX 



collaborations 16, 17 



In the present calculation we have restricted to lowest-order QED calculations paying a 
special attention to realistic form factors and absorption effects and totally ignored higher- 
order corrections. How important are the QED higher-order correction was demonstrated 
recently in Refs. 0, [l3| • While Jentschura and Serbo [13] argue that the higher-order cor- 
rections are rather small, Baltz [t^ finds a huge reduction of the integrated cross section of 
the order of 20%. Cleary the discrepancy should be clarified in the future. It would be also 
very interesting to calculate the higher-order corrections for differential distributions which 
will be measured at LHC ^. The latter calculations seem to us rather difficult technically. 

In the moment it seems precocious to answer the question whether the processes dis- 
cussed here could be used as a luminosity monitor for heavy ion collisions at LHC. In our 
opinion, first these processes should be measured and compared to theoretical calculations. 
In addition, the infiuence of the absorption effects and multiphoton processes on differential 
distributions should be studied in more detail. 
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